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ABSTRACT 


The calculation of Viscous flows with coupling methods 
is surveyed. The approximation levels and concepts are 
first outlined, as well as the generalized formulations of 
the viscous displacement over the inviscid flow. Then, 
the strongly interacting methods that are based on thin 
viscous layers approximations are discussed. 

In this way, a matching formulation of the viscous 
flow, calculated as a difference with the inviscid overlay- 
ing flow, is suggested/ in order to restore approximately 
the normal pressure gradient inside of the layers, as well 
as to remove the supercritical behaviors, in the Crocco- 
Lees sense. This analysis maintains simple viscous inte- 
gral equations. A review is then given of the main numeri- 
cal techniques presently available for the coupling problem. 
The global state of the art and possible extensions are 
looked at through the viscous methods for airfoils. New 
results are presented for trailing-edge separation, and 
an approximate method is suggested to capture numerically 
the viscous interaction under the shock-waves. 
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COMPUTATION OF FLOWS INCLUDING STRONG VISCOUS INTERACTIONS 

WITH COUPLING METHODS 

J.C. LeBalleur 


I . INTRODUCTION 


High Reynolds numbers flows are characterized by the /l* 

existence of very rapidly changing localized phenomena that 
closely associate the effects of viscosity and turbulence, 
generally contained in the thin layers, and their effect on 
the structure of non-dissipative flows. A typical example 
of such a condition is the shock-boundary layer interaction. 

These concentrated arid small-scale phenomena may con- 
dition the macroscopic flow completely, owing to the non- 
linearity of the equations which govern them and which gen- 
erate complex flow conditions, called strong viscous inter- 
actions, in contrast to flows with weak interactions which 
are fundamentally different and lead to solutions close to 
conventional approximations of inviscid flows. 

Notwithstanding the considerable progress made in adap- 
ting numerical methods to the resolution of complete Navier- 
Stokes equations, which in principle are able to process 
strong interacting flows, a purely mathematical approach 
to the aerodynamic problem remains difficult to the extent 
where it remains highly dependent upon the accuracy of ; 
the. numerical techniques, the preparation of large enough 
turbulence models, the necessary recourse to very fine 
meshworks requiring the use of very power, and still' 'un- 
available, computer means. 

Accordingly, a more physical approach is present in 
all practical numerical aerodynamics. It appears, among 
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other places, when introducing equations for boundary 
conditions or for simplified local "models", when pre- 
paring composite resolution methods based on the coup- 
ling of several distinct numerical methods, when it is 
necessary to use meshworks adpated to local flow condi- 
tions. It is clear that the origin of such methods is 
intricately interrelated with the multiple-scale struc- 
ture of the flow as well as with the existence of approx- 
imated equations and numerical methods adapted to each 
scale. The adaptation of these methods to the physics 
of the phenomena has often proven to be a source of ef- 
ficiency, particularly for relatively modest levels of 
viscous approximation. A typical_ example is given by 
the difficult calculation of transonic airfoils for which 
methods based on inviscid flow and interacting boundary 
layer analyses exhibit even better performances than the 
Navier-Stokes type analyses, both in terms of computer 
costs and their accuracy. 

On the long-term, we may conclude that the adapta- 
tion of coupling methods to the specific physical aspects 
of high Reynolds numbers should lead, in conjunction with 
the development of direct resolution methods, to greater 
economy for equal performances, or to the processing of 
more complex problems for a given computer cost (three- /2 

dimensional flows, for example) , or to an indirect devel- 
opment of new numerical techniques for resolving Navier- 
Stokes equations. 

Parts 2 and 3 of this report will succinctly examine 
the approximation levels and the general concepts which 
seem to be outlined from current methods of computing 
flows with strong viscous interactions. Parts 4, 5 will 
discuss methods based on simple viscous approximations, 
obtained from Prandtl or -thin layer hypotheses. After 
discussing their effects on a computer model by coupling 
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with the inviscid flow, we will briefly survey the numer- 
ical methods currently capable of resolving a coupling of 
this type, which is a basic i problem whose study is still 
very recent and limited. Finally, we will take a lo6k at the 
capabilities and future prospects for the coupling methods 
of this category by examining an example offering a synthesis 
of transonic airfoils where various interacting viscous phen- 
omena interfere. 

2. NUMERICAL APPROXIMATION OF VISCOUS FLOWS: GENERAL 
LEVELS OF ANALYSIS AND COUPLING METHODS 

It is not within the scope of this brief analysis to 
review the various methods used to calculate viscous flows, 
as our objective here is simply to offer a qualitative out- 
line of the main 'approximation methods, the possible approx- 
imation levels, and their application limits. 

As shown in figure 1, a classification into four lev- 
els seems conceivable^ level I combining the various meth- 
ods based on a direct resolution, while levels II, III, IV 
correspond to the main indirect approaches or coupling meth- 
ods. Although relatively strong approximations may be used 
for each level, which prohibits a rigorous structural organ- 
ization, the general nature of the most complete methods of 
each category increases from level IV to level I. 

2.1. Methods of Direct Resolution 

Level I corresponds to the methods of direct resolution 
and includes in particular the complete resolution of the 
Navier-Stokes equations, as well as the resolutions limited 
to large-scale turbulent structures, with modelization of 
fine structures. From a standpoint leaning toward methods 
currently usable in applications, level I corresponds to an 
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approach which may be qualified as a general approach [1 to 
4J to the extent that it adopts the same numerical method 
and the same system of equations which are uniformly valid 
in the entire flow field. We are referring to averaged 
Navier-Stokes equations, associated with a turbulence model, 
or Navier-Stokes equations truncated by thin layer approxi- 
mations [5,6]. 

2.2. Direct Methods With Equation Commutation 

The first extension of the overall approach consists 
of retaining a single numerical method, but by commutating 
the resolved equations over the complete or simplified sys- 
tems 11 ], figure 2, as a function of the calculation regions. 
This approach, similar in its intentions to the methods of 
coupling by regions, cannot be regarded identical to it. It 
is actually clear that the equation commutation boundaries 
thus employed do not in any way materialize distinct regions 
of numerical calculation, linked by their boundary conditions 
alone; on the contrary* a numerical . continuity of the deri- 
vatives at the commutation boundaries is implicitly postu- 
lated. The approach of commuting equations by region es- 
sentially consists of simplifying the calculations by assum- 
ing a priori that part of the equation terms are discarded 
in certain regions, in accordance with the precision of the 
numerical resolution technique. This simplification may be 
completed by optimizing the geometry of the commutation 
boundaries . 

2.3. Str.ohg Coupling Methods With Division Into Computation 
Regions 

The equation commutation approach may generate new 
mathematical or numerical problems compared to an overall 
approach with complete Navier-Stokes equations, in the event 
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of a large truncation of these. For example, they will 
be associated with order changes in equations with par- 
tial derivatives, changes in the domain of dependence and 
in the required boundary conditions. 

This leads us to level II and to the coupling approach 
by separating the computation regions (figure 3) in which 
the mathematical and numerical problems are treated separ- 
ately, while the overall problem is restored by means of 
a rigourous coupling, called strong coupling, between the 
boundary conditions of the various sub-regions along their 
common boundaries. The equations used in the various sub- 
regions are generally either complete Navier-Stokes equa- 
tions, or Prandtl or :Euler\ 'equations truncated in thin 
layer approximations [8, 9, 10, 11]. 

Changes in the order and nature of equations with 
partial derivatives with computation sub-regions forces us 
to modulate the number of coupling conditions along the com- 
mon or connecting boundaries, and to make selections on the 
basis of the physics of the phenomena as well as on the 
nature of the local equations. The coupling of aerodynamic 
parameters to a connecting boundary no longer leads in this 
particular case to a linking of their normal derivatives to 
this boundary, the lack of residual connection being repre- 
sentative of the precision of the coupling model, for a 
determined division of regions and equations. 

Note that the transition from level I to level II 
creates a new numerical problem. This consists of making 
a strong coupling of boundary layers associated with various 
sub-problems, using iterative or implicit methods in steady 
or unsteady conditions. It is clear that this major pro- 
blem could not be summarized by establishing a computation 
flow chart, and in aerodynamics it has been the subject 
matter of a limited number of studies essentially motivated 
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by linking the Euler and Prandtl equations (see paragraph /3 

5) . 

2.4. Methods Of Strong Coupling With Division Into A 
System Of Equations 

The technique described above of the division into 
domains supposes that in each of them the equations re- 
solved themselves describe the local behavior of the flow. 

In the case of modest approximation levels, like that of 
the Prandtl equations, such a condition may lead to large 
limitations. It should be remarked that such a connecting 
technique' is to be avoided in the analytical methods of 
connected asymptotic expansions, in which the calculation 
domains associated with the various scales' become super- 
posed in the same physical domain. 

We have found [3, 12 to 15] that the introduction of 
overlapping calculation domains , (figure 4) in the formula- 
tion of coupled composite methods for calculating viscous 
flows at high Reynolds number, could lead to considerable 
progress, physically related to the possibility of proces- 
sesing the viscous flow in terms of differences with the 
fictive inviscid flow which it is the closest to. This 
analysis constitues a generalization of the boundary layer 
methods, whether conventional or with strong coupling along 
a wall, which we will assume to be plane in order to simpli- 
fy the presentation, if w and v are the components for the 
> velocity, p the density and p the pressure, if u 3 v 3 p, p 
are their homologues for the overlaid inviscid flow, the 
equations for the viscous flow, notated j(T) + (?) will be dis- 
sociated into a system of inviscid flow equations, notated 
j(j) and into a system of matching viscous equations notated 
Kz).~ The unsteady dynamic equations therefore become: 


6 



(Z.l) 



dpu . dpir n 
dx dy 




; © ' 

dpu 1 dpJir dp 

dx dy dx 




i 

dputr dpi ) 1 dp. 

dx dtj djf 




© < 

dpa-pu. dptr-pir q 

dx 1 dy 

dpMd.pu . 1 dpMO- pMtr 

dp.p 

i 

j 

i 

dZxq 

dx ' dp 

c>x 

dx 


\ 

dpMV . putr dptr pn ~ 

dp.p 




dx 1 dtf 


* 

dx 


o x3 d U i C Xv represent the stress tensor components. This 
matching formulation of the viscous problem leads to a 
very similar resolution as in the method of coupling by 
domains, owing to the strong coupling that must be achiev- 
ed between and the boundary conditions so that: 

0 = lim [f-f] with f = {ujV 3 p 3 p} (2.3) 

We may still notice that the importance of selecting 
the external boundary of the viscous layers y = 6 (x ) , 
which is relatively arbitrary in the method of coupling by 
linking the domains, is eliminated. 

An advantage of the matching formulation (2.2) appears 
clearly by approximating with intermediary assumptions 
which fall between Prandtl's assumptions and those of Navier- 
Stokes equations for thin layers. We therefore have: 
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which, according to (2.3), gives p (x>y) =p(x*y). 


This approximation of the pressure is much less Jr.es trie- /4 

tive than that deduced from the Prandtl equations in a 
coupling by linking the domains: 


M-*0 f,(x)m f. (X,S) 

Oy ’ 


(2.5) 


while the resolution of the parabolic system (2.4) with 
«, u., p, p,. given by (2.1) is as complex as the resolution 
of simple Prandtl equations. 


In the event of less Extensive simplifications than 
(.2.4), or of an application of the Navier-Stokes equations 
(2.2 )l, the selection of the matahing formulation becomes 
intricately related to the selection and conditioning of 
the numerical techniques applied to (2.1) and (2.2) as 
well as to their coupling. Finally, we will note that at 
approximation level II in figure 1, the division made into 
a system of viscous and non-viscous equations implies an 
analysis of a strong interacting coupling, which corres- 
ponds to a strict resolution of (2.1) (2.2) (2.3), and 

which accounts for the influence exerted by the downstream 
flow on the upstream flow. Approximations of this type, 
although relatively simplified, have been studied recently 
116, 17, 18J . Conversely, other approaches maintain ap- 
proximations of systems of parabolic equations , resolved 
directly from upstream to downstream, which makes them 
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comparable to boundary layer methods for weak interactions 
of level IV 119, 20]. 

2.5. Weak Coupling Methods 

Approximation levels I and II correspond to a model- 
ization of the viscous flow capable of describing flows 
of strong viscous interactions, at least qualitatively, 
owing to the viscous equations resolved or to the strong 
coupling effect at the inviscid flow, in the case of the 
simplest viscous equations. We should still consider that 
this capacity of processing regions of strong interaction 
implies the use of suitably fine interconnections in order 
to be effective, for any equation resolved. In turbulent 
flows, even for the simplest systems of equations, the 
resolution scales should not be coarser than the local 
thickness of the viscous layers involved. 

Levels III and IV correspond to simplified analyses 
whose application generality is smaller, but which remains 
the support of most current applications. The most com- 
monly used are the boundary layer methods of level IV, used 
in conjunction with conventional inviscid flow calculations 
resolved in uncoupled mode. The absence of a coupled reso- 
lution comes from the concept of a weak interaction or coup 
ling, related to the asymptotic second order theory of the 
boundary layer in laminar cases [21] , and to theories of 
the same nature in turbulent cases [22,23]. 

Figure 5 shows the scheme of the' simplification 
made at the coupling level. In a first approximation, the 
field> of inviscid flows is obtained independently of any 
boundary layer. It determines a simplified pressure field 
for the boundary layer which makes it possible to calculate 
the boundary layer and its displacement thickness in the 



first approximation. This makes it possible to determine 
the final approximation of the inviscid flow, as well as 
that of the boundary layer pressure field. A second order 
approximation may then be calculated for the boundary lay- 
er, which then always has a normal predetermined pressure 
gradient, but which is not zero in the case of a wall cam- 
ber. The coupling is "weak" provided that the viscous lay- 
er is completely conditioned for a given approximation order 
by the inviscid flow upon which no mutual influence is ex- 
erted . 

In this analysis, the notion of iteration on the vis- 
cous displacement thickness is absent and is replaced by 
a direct and single correction. Any shortcoming in such 
a single correction precisely marks the appearance of phen- 
omena of strong viscous interactions, related to the separ- 
ations, shock waves, training edge, as well as to the ne- 
cessity for multiple layered models if we still want to 
use asymptotic theories [23 to 31] without being able to 
eliminate the appearance of a strong coupling problem. 

A characteristic property of a weak asymptotic coup- 
ling is that it leads to the resolution of a parabolic 
system, for the viscous problem, in which we do not know 
the effect exerted by the downstream flow on the upstream 
flow. Such- an influence can be determined only through 
the boundary conditions set by the inviscid flow, and 
which is therefore totally eliminated for a supersonic 
external flow, or even in the approximation of a one- 
dimensional subsonic inviscid flow. 

The corresponding simplification made of the resolu- 
tion can therefore be carried out upstream to downstream, 
and has been the subject of transpositions in methods 
abandoning the simple Prandtl equations for less restric- 
tive equations, qualified at level Iv in figurel of the 
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parabolic approximations of the Navier-Stokes equations. 

In this case, the Navier-Stokes equations were first 
subjected to thin layer approximations, which consisted 
of eliminating the viscous diffusion terms according to 
Ox, if Ox represents the direction of the main flow. This 
simplification of levels I or II is not self-sufficient in 
the general case where the system of equations is complet- 
ely parabolized in the direction of positive x's, particu- 
lary in the presence of subsonic flow regions [32, 33] . 

We may therefore speak of systems of elliptical equations 
because of the pressure field, provided that the systems 
remain parabolic either when the pressure is fixed, or 
when its mean longitudinal variation along x is fixed, 
while its secondary variation in plane yOz is calculated 
134] , or when the longitudinal pressure variation is also 
calculated upstream to downstream, but at the cost of a 
quasi-one-dimensional approximation and of a pressure 
change in the momentum equation according to Oa: [35] . 

In all of the cases, the switch from such a parabol- 
ized weak coupling problem to an elliptical strong coup- 
ling problem of level I or II, implies the addition of an 
iteration on the pressure, capable of eliminating at con- 
vergence the various approximations used for each parabol- 
ized iteration [16, 17, 34, 36]. 

2.6. Intermediary Methods Between Weak And Strong Couplings /5 

Figure 1, of level III, shows examples of this type 
of analysis. The first may be given through an progressive 
extension of the methods of computing airfoils. Starting 
with the obvious short-coming of a weak coupling approach, 
an extension toward the investigation of a strong coupling 
was progressively developed, by means of iterations on the 
viscous displacement effect. Although the approximation 
levels obtained differed considerably between the different 
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methods, for example 115, 29, 37 to 40], either because 
of the use of smoothings or relatively arbitrary modeli- 
zations of the displacement thickness, or due to numeri- 
cally or analytically incomplete coupling techniques in 
regard to separations, trailing edges, boundary layer/shock 
wave interactions or calculation meshings, it is clear that 
the overall variation should lead to strong coupling meth- 
ods of level II. 

A second category of intermediary methods between 
weak and strong couplings is given by forming a model of 
inviscid flow separations. A first example is given here 
by calculating inviscid flows containing vortex sheets 
issuing from predetermined separation lines [41] . A second 
example is that of forming models 6f large separations by 
calculating free isobar boundaries in an inviscid flow. 
These models are especially used in stall problems [42] . 

3. BOUNDARY CONDITIONS OF INVISCID FLOWS - DISPLACEMENT 


The methods of calculating inviscid flows is based 
on the slip condition expressing wall tightness, arid on 
the Joukowski condition which when used seems to recall 
alone that the objective is to form a model of a viscous 
flow. The calculation also brings to light, in regard to 
weak solutions, the slip sheets through which the normal 
velocity and pressure are continuous. 

If the physical problem is approached by analyzing 
a viscous flow, at the scale of a continous medium, such 
slip conditions become impossible for any Reynolds number 
and it is the asymptotic theory of the boundary layer 
which tells us that the external non-viscous flow to be 
considered is subjected to these slip conditions, in a 
first approximation. These slip conditions therefore 
acquire an origin related to the viscosity. At the 


second approximation order, a normal non-zero velocity on 
the walls, as well as a normal velocity and pressure dis- 
continuity on the slip sheets, express the viscous displace- 
ment effect on the residual inviscid flow region. If in- 
dices 1 and 2 designate the successive orders of approxima- 
tion, we have for example: 



1 ■ a ''/ 


& < Oleyno/Ji 

! r M [ e<“J(x,oj -* a tJjlV’lfx.o) - 1 J d $ 

k r *ho) 


( 3 . 1 ) 


3 . 1 . strong Coupling Along Adjoining Regions 

The transposition of these conditions to coupling 
problems at finite Reynolds number was initiated by Crocco 
and Lees 143] to calculate supersonic separations in a 
strong coupling approach along adjoining regions (figure 6) 
using the Prandtl equations for 0 <y<S(x)> S being the phys- 
ical thickness of the boundary layer. In the most usual 
case where only one boundary condition of viscous origin 
in y = S(x) is needed for determining the inviscid flow, 
the relationship used by Crocco and Lees is general, since 
it is based only on the integration of only one equation 
of continuity in y. If K(x) is the wall camber: 
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ft*) [ e“J(x,j) * / (M(x.s) ~[e*](x.y)j % 


( 3 . 2 ) 


This displacement formulation 8*(x) is exempt from approx- 
imation other than that associated with the calculation of 
[py] (x'j y ) . 

Other methods of introducing the coupling condition (3.2) 
in the calculation of inviscid flows are frequently used. It 
may seem more convenient to consider that there is either no 
flow through the boundary of the inviscid flow, and we are 
then speaking of a displacement concept, or that this bound- 
ary is simply based on the wall. In the two cases, a cer- 
tain equivalence with (3.2) may be obtained by considering 
an approximated analytical shift of (3.2) from y = 8(x) in 
y = y(x). The equation of non-viscous continuity gives: 


ll-ky, ■ (i-ko - (y - S H)r %-fcr, * - (3.3) 

which leads to a slip condition if y = 6*, an injection /6 

condition if y = 0 and if the corresponding expressions 
for 5 * (x) : 


ftr } dS' 

hrJ/xj') = ' -• 

( 3 . 4 ) 

Ct*]fx.6) = 4r s '](*.o) * 

( 3 . 5 ) 

s MhJ ( x.y ) -]' {[e*] (x . Y) -[?Z](x.y) )d, 

( 3 . 6 ) 


The displacement formulations (3.4) (3.5) (3.6) are still 

approached, as well as the equivalence with (3.2). The 
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asymptotic validity of this equivalence when <?— « is not 
automatic where the expansion order (3.3) depends on the 
order of dp]x/dx, a quantity often not bound when 
in the asymptotic theories with multiple layers. The 
formulations (3.4) (3.5) may still be outlined from any 
approximation by selecting new definitions for 8* (x). 


3.2. Stong Coupling Along Overlapping Regions 


Such a change in the definition of 8*(x) is easily 
done by introducing (figure 7) overlapping calculation 
regions over which the viscous equations ’.are segmented ’ 
as described in paragraph 2.4. In this case, the simple 
integration of the continuity equation (2.2) in y gives, 
since 

o - e v ] • 

M(x.o) - 

f(z) [e*] (x.0) » J t {[H(z,y) } d t 


(3.7) 


\ 0, lim [f 

! a -*a L 

I ' 


This generalization of the definition of 8(x) takes into 
account the normal pressure gradients inside the viscous 
layer, through variations in y of the viscous term cnT and of 
the non-viscous term py; these may possibly be rapid, even 
discontinuous, unlike boundary layer analyses. A schemat- 
ic illustration is given in figure 8, for the case where 
(pT/y) = $(Xjy) is assumed to be a continuous function. 

A formulation of 8*(x) adapated to the displacement 
concept may likewise be defined: 
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J [e* ](*.y) d ! - Jj ?“](*,?) d D 


( 3 . 6 ) 


All things considered, use of a displacement thickness to 
formulate a coupling problem is in no : way the same as us-t 
ing boundary layer approximations. 

4. SIMPLE APPROXIMATIONS OF THE VISCOUS FLOW - INFLUENCE 
OF A STRONG COUPLING 


For both historical and practical reasons, current 
coupling methods have been mainly developed on the basis 
of boundary layer assumptions. We will confine ourselves 
in the remainder of this report to this class of methods, 
using either Prandtl -equations, or a similar extension of 
thin layer equations, described in paragraph 2.4, in which 
the pressure p(x,y) is identified with that of the non-vis- 
cous calculation p(x 3 y). 

4.1. Separation and Reverse Flow Regions 

The boundary layer problem, subjected to the Prandtl 
equations and decoupled from the inviscid flow, must be 
formed by an external boundary layer condition, for which 
the pressure p(x>&) has been used for a long time in the 
same way as in the asymptotic theories of weak interactions. 
We thus know the singular behavior of Goldstein [44] for 
"square roots" offered by solutions near a position of zero 
friction in two-dimensional steady flows, a behavior also 
offered by integral boundary layer methods [12] . In un- 
steady flows 1 4 5 j or three-dimensional flows [46], the 
integral methods demonstrate a different type of singu- 
larities which are formed of discontinuities of integral 
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thicknesses associated with the weak solutions of equa- 
tion systems, which are no longer initially localized in 
reverse flow regions, and which might explain certain num- 
erical anomalies found in the solution of the Prandtl un- 
steady equations with reverse flows and a set pressure 
field. 

However, we were able to show that these singular /7 

behaviors should not be confused with a validity limit of 
the Prandtl equations, or with an abrupt separation from 
the viscous layer and from the wall, but that they result 
from the^pressure selected for the external boundary con- 
dition. When the pressure becomes a calculation unknown, 
either by formulating an inverse problem subjected to an 
external condition, or by coupling with the inviscid flow, 
consistent and realistic solutions were obtained for small 
separations [12, 47 to 55]. The same behavior Was bbserved 
for unsteady conditions [45] . YetJ.t still clear that that 
the Prandtl problem for strong couplings with the external 
inviscid flow no longer constitutes a parabolic system 
free of downstream effects, even without the reverse flows, 
owing to the boundary condition, unknown a priori, formed 
by the external pressure distribution with coupling 'ef- 
fect [figure 6). 

4.2. Downstream Influence On The Upstream In Supersonic 
And Transonic Flows 

The essential question raised by the matching calcu- 
lation of the Prandtl and Euler equations is the following: 
is thh '. elliptical nature of the unsteady Navier-StOkes 
equations also~ found in the coupling model when the 
inviscid flow is locally supersonic so that the down- 
stream influence on the upstream cannot have an exclus- 
ively non-visous origin? Since Crocco-Lees' first work, a 
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a long investigation was conducted [12, 14, 56 to 59] and 
the answer found was positive for "subcritical" viscous 
layers, and negative for "supercritical layers.'! 

In the subcritical case, the coupling of hyperbolic 
non-viscous equations with parabolic viscous equations 
leads 'to a problem of initial conditions which is improperly 
presented, to the extent where an exponential amplification 
of any initial disturbance occurs on a short scale, corres- 
ponding to the boundary layer thickness, approximately, 
according to a branching solutions process. The aspect of 
an improperly presented problem is eliminated if one of the 
boundary conditions is transferred from upstream to down- 
stream, which reveals the elliptical nature of the model. 
Furthermore, the branching process generally makes it pos- 
sible to introduce such a downstream boundary condition as 
a purely optional constraint which is the case in the coup- 
ling 'model for the boundary conditions currently applied 
to the downstream in Navier-Stokes type calculations. The 
occurrence of a down tr earn to upstream influence due to the 
coupling effect also finds a theoretical support in strong 
interaction asymptotic theories, such as the triple-deck 
laminar model, for example. We were also able to verify 
numerically that the supersonic triple-deck equations are 
also asymptotic solutions of the coupling model, when the 
Reynolds number moves toward infinity [27], figure 9. 

In the supercritical case, the elliptical nature of 
the viscous equations is only partially found in CtocCo- 
Lees coupling model, by means of weak solutions whose 
pressure jumps suddenly, transform the supercritical lay- 
ers into subcritical layers, each time a downstream in- 
fluence on the upstream is locally indispensable. The 
supercritical behavior is '.therefore found agains near 
the downstream after passing a "critical point," where 
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the stability required of the solution plays the role of a sup- 
plementary downstream condition for the subcritical region [59] . 

A simple interpretation of the behavior patterns of the model in 
the vicinity of jumps and critical points may be directly deduced 
from the local behavioral patterns of the branching solutions, as 
well as the stabilizing boundary conditions needed, see figures 
10 to 11. 

The preceding gaps in the Crocco-Lees model are not related 
to the approximations used to solve the Prandtel equations if an 
integral method is used, as was often believed to be the case, 
but are actually due to the simplified pressure field p (x) = 
p(x,S) of the viscous layer. In an integral method the rela- 
tionship (3.2) leads to a quasi-linear relationship between the 

.f * . — 

viscous deflection ; &-(%.) and the pressure gradient in y =-S(x): 

\D, 0 » D z S*±. , D s (4.1) 

vanishes to zero for separations or for reattachments where the 
pressure gradient cannot be any value, but must be S*(dp/dx) = 
-(C^A^) to assure the regularity of the solutions. reduces 

to zero at the critical points, where the regularity forces 
9. /Di ) . The symmetry of the Goldstein and Crocco-Lees 

singularities may be noted for the coupling model for adjoining 
regions. The relationship (4.1) or simply 

!Q . 3 4L + C B ~ -7T (4.2) 

j dz Vi 

is a determinant for the branching solutions, which are stable 
according to x if B < 0, and unstable if B > 0. For an attach- 
ed layer, the unstable case corresponds to a subcritical be- 
havior and to an elliptical type problem. Note, however, that 
in the case of a separation, D ^ changes its sign, and therefore 
so does B, and this leads to table branching solutions in the 
separated region and the downstream effectuon the upstream^has the 
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result of selecting stable solutions at the reattachment 
point, figure 12. In the case of a local solution of the 
Prandtl equations, we are left [57] with a relationship 
analogous with (4.2), demonstrating the importance of the 
local Mach number M in supercritical behaviors: 



(4.3) 


The limit e -*■ 0 often being singular, the behavior of the 
branching solutions is generally not deduced directly from 

(4.3) , see [3, 12, 15]. Nevertheless, it follows from 

(4.3) that the supercritical behaviors (B < 0) can be de- 
rived only from the presence of supersonic viscous regions, 
and that they are highly dependent upon the selection of 
the boundary 6 in the matching model. Any subcritical 
coupling may thus be converted into a supercritical coup- 
ling, by selecting a boundary y = S(x) farther away from 
the wall [3] . 


The subcritical or supercritical nature of Crocco-Lees 
therefore iis not an intrinsic property of the viscous layer, 
but of the ad joining coupling model* [3, 15]. Its physical /8 

interpretation should therefore call on the external in- 
viscid flow. We may thus note that the supercritical be- 
havior of a supersonic separation at increasing Mach number 
appears only when the focalization shock wave approaches 
until it reaches the external boundary of the viscous layer, 
with any penetration inside the viscous layer being excluded 
by the calculation model and being simulated by a pressure 
jump, figure 13. 
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4.3. Approximated Treatment Of the Normal Pressure Gradient. 

Matching Formulation 

The shortcomings just described of the simple Crocco- 
Lees model for processing flows with strong viscous inter- 
action in supersonic states, as well as in transonic states 
(M ±1.201, are the result of a too simplified pressure 
field pCx) = pCxjSJ. However, we have found [3, 12, 13, 15] 
that it is quite easy to discard these shortcomings by 
establishing a matching formulation of the viscous calcula- 
tion, corresponding to the viscous calculation regions and 
the overlapping non-viscous regions, figure 7, as described 
in paragraph 2.4., and if the approximation over the pressure 
becomes p (x 3 y) = p(x 3 y). 


By retaining only Prandtl's viscous term in the system 
(.2.2), we obtain by integrating in y the equation of motion 
along x and its first moment relative to u: 


M ft**'**- ^ d 0' M(x - o) " Z(x ’ 0) " 4 

f[e** *j[ d ? 


(4.4) 

(4.5) 


while the equation of motion in [x 3 SCxJ] gives the 
entrainment equation: 


with \ dS_ _ f ^ « £ (4.6) 

' j * ’ 1 ? f * £ '{ Lo 

The approximation p(x 3 y) = p(x 3 y) also discards the 
last term in (4.4) and in (4.5). If we define the general- 
ized integral thicknesses for motion and kinetic energy in 
an analogous mode to that adopted for the displacement 
thickness in (3.7): 
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(4.7) 


[s%) t b (X )][^] m - £[t**‘e* z ](».}) d 2 

[ 3%) * &*(x)] [ ’e**](x. o) » £ [p* - e*’J ( x ,tf ) d y 
i*i(x) A( Zi o) « J' [ji-Z ]( x ^ dy 


We obtain integral equations very close to those derived 
with the Prandtl equations: 


l: 

h 


d& 3* * 

dx * * 


ZB du 9 df J Cf 

Jz + P dz ~ 2 


JL M. + z IzJL L o) = (p * 4> 

dx * pa* dx *. dx 


(4.8) 

(4.9) 


Note, however, that the kinetic energy equation (4.9) 
requires a modelization of a term for the normal pressure 
gradient <pp , related to the variation of 9 p/dx in y in (4.5), 
whereas the entrainment equation (4.6) does not require it. 

A certain approximation is still present, however, if the 
entrainment equation (4.6) is transformed so that only the 
following terms occur: p (x,0) 3 u(x,o), v(x,o): 


/ f dS _ d!" 3.5 * dpu J £ 

I [ dx dx pa dx J( X '°J 


(4.10) 


Besides the benefit of accounting for an approximated 
normal pressure gradient, the matching formulation (3.7) 

C4 . 7 to 4.10) leads us to replace the integral B in (4.3) 
with an integral B f £12] for the analysis of supercritical 
behaviors: 

| u l\ i dn * nonviscous 

j £ -° 4 3 T- ** fi(y) - viscous (4.11) 


The first improvement, commdn to (3.7) (4.7 to 4.11) 




22 



is to totally discard the relatively arbitary role of 
the adjoining boundary y = 8(x). Secondly, the term 
(1_M Z ) in (4.3), generator of supercritical behaviors 
when M>1 , is replaced by (A / 2 — Af 2 ) , a term which is most /9 

likely always positive. We were thus able to conclude 
[12] that the coupling along overlapping regions always 
leads to subcritical behaviors, where the downstream in- 
fuelnce on the upstream is fully observed, and where the 
pressure distributions on the walls p(x 3 o) are always con- 
tinuous. Consequently, the possible shock waves always 
generate by focalization within the non-viscous calculation 
region, which may even occur within the viscous layers, fig- 
ure 14. in these intense cases of normal pressure gradient, 
a second order may be applied a posteriori to the approxima- 
tion p (Xj y) = p(x 3 y), , in the form: 

! £ . k (*) 

L dj 

K(x) is a mean camber of the stream lines, which may be 
indistinguishable from that of the walls or mean lines of 
waves only in regions of weak interactions. At the wall, or 
at the center of a wake, we obtain: 

pCx 3 o) - pCx 3 o) = K(x) = Q(x)] (4.13] 


[eu 1 - e*](x.y) 


(4.12) 


Finally, it should be pointed out that the super- 
critical behaviors, eliminated in the preceding overlap- 
ping coupling should still be possible in the coupling 
by adding a displacement thickness, where relationship 
C4.ll) becomes: 





M 1 . M l _ 
Vf. M l M 1 


d f 


(4.14) 


an integral in which the first term may always be negative. 
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threshold of transition to supercritical behavior 
therefore being essentially M = 2 for a flat plate tur- 
bulent layer . 

5. NUMERICAL METHODS FOR STRONG COUPLINGS 

The simplifications brought to the resolution of the 
viscous fluid, at level II in figure 1, by segmentation 
of the calculation regions or of the equations used, have 
another role of artificially generating a strong coupling 
numerical problem. In the case of a coupling by adjoining 
regions, it is necessary to assure a rigorous compatibility 
of the boundary conditions of the various sub-problems on 
their common boundaries. Since the nonlinearity of the 
viscous equations in practice makes an iterative solution 
inevitable, it may be assumed that the problem remains 
virtually the same in the coupling by segmentation of the 
equations, with overlapping calculation regions. 

5.1. Nature Of The Problem 

This strong coupling of the boundary conditions over 
the interfaces of the various constituent calculations re- 
quires a less rigorous numerical resolution than that used 
in each subproblem, with the risk of degenerating the 
strong coupling treaetment into a weak coupling approxima- 
tion. This condition is reinforced by the high sensiti- 
vity of the inviscid flow calculations to the accuracy of 
the numerical treatment of the boundary conditions. It 
is even more sensitive if we consider that in thin layer 
viscous problems, the coupling effect alone can change 
the mathematical nature of the behavior of the various 
subproblems, via branching solutions reviewed in paragraph 
4.2 and lead to differences of about 1 between the solu- 
tions subjected to a rigour coupling and those derived 
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with an appromixated coupling. 


If the finite differences methods are used, for 
example, the coupling must be made in each discretiza- 
tion node of the interface with the consistent numeri- 
cal schemes, compable to those used in each region. The 
interpolation over the interface is conceivable in the 
case of discontinuous meshes, provided that it is consis- 
tent with the numerical scheme order, and with the regions 
of mathematical influence. Conversely, the use of smooth- 
ings or filtering is to be discarded, unless the technique 
may be reduced to the insertion of a numerical dissipation 
which is evanescent when the discretization step leads to 
zero . 


If we consider the solution of evolving problems, the 
same problem with a strong coupling may be overcome at 
each time step to avoid a viscous resolution inconsistent 
in time or a weak interaction approach. In the case of 
iterative relaxation or pseudo-unsteady methods, in which 
only the final solution matters, the stability in time may 
be disregarded. Conversely, coupling errors at each iter- 
ation should not augment and a specific stability study 
should be considered. The specificity and difficulty of 
such a study is due to the fact that although it is nec- 
essary to analyze the magnitudes bbrdefing over the coup- 
ling interface, error increases is managed by an operator 
calling on the resolution of the subproblems adjacent to 
this interface, which defines a mathematical problem over • 
a space whose dimension is larger than that of the inter- 
face where the coupling errors are studied . 

Yet, it is clear that, apart, from unusual circumstan- 
ces, the definition of an iterative coupling method may be 
reduced to a simple definition of a calculation flow chart 



in the case of a strong coupling. In the first place, 
it is necessary for the coupling to select consistent 
numerical schemes, at least in terms of space, observing 
not only the dependency domairis of the ^equations of 
each subregion calculated separately, but also the real 
dependency domain of the coupled mathematical problem, 
specifically accounting for the possible branching solu- 
tions. For example, the consistency will not be obtained 
for the coupling of the Euler and Prandtl equations along 
a locally supersonic boundary, with attached buscritical 
boundary layer, if the schemes with differences totally 
decentered upstream are used and if, at least for quanti- 
fying the coupling relationships, one of them does not 
escape this choice. On the contrary, we have shown in 
this case, the advantage of a decentering toward downstream 
[13] . 


In the second place, the generality of an iterative /10 

coupling method implies that the complex stability problem 
is mastered, for the purpose of adapting the relaxation 
techniques, i.e. modulating the stabilizing relaxation 
coefficients as a function of time intervals and local 
spaces, as well as of the aerodynamic magnitudes. This 
constraint is similar, to the necessity of calculating, 
in an explicit unsteady numerical method, the local maximum 
time interval compatible with the stability of the calcula- 
tion as a function of the local space interval. 

The elimination of this constraint has for a long 
time led to unstable iterations for boundary layer calcu- 
lations over airfoils and to a relatively, arbitrary use of 
smoothing techniques. These instabilities, which are purely 
numerical, are influenced by the weak or strong aspect of 
the local viscous coupling interaction, but they are not 
always absent if the discretization interval is small 
enough. They particularly appear for flat plate boundary 


26 



layers [60] . Figure 15 gives an illustration relative 
to the inviscid flow - unsteady boundary layer coupling 
[61] near a uniform flow over a flat plate. It shows 
that, all other things being equal, a slight reduction 
in the space interval or increase in the time interval 
is sufficient for switching from a stable 1 qoupling iter- 
ation to an unstable iteration. 

Finally, a reduction' of the numerical instabilities 
of explicit type simple ^iterative coupling techniques, 
may obviously be investigated in the development of 
methods which process the coupling in a more implicit 
manner, for each iteration, or based on more general 
nonlinear numerical techniques. 

5.2. Behavioral Law of The Steady Boundary Layer 

Most of the viscous - nonviscous numerical coupling 
methods have a slightly implicit nature, for reasons of 
simplicity, application generality and interchangeability 
of the viscous and nonviscous modules. Recourse to vis- 
cous modules using Navier-Stokes equations, which are com 
plete or in the approximation of thin layers, remains 
rare in the coupling methods [8 to 11, 32 to 34] for the 
part specific to the coupling. Conversely, more accurate 
methods are outlined in the case of simple viscous approx 
imations as described in paragraph 4. 

In the case of integral boundary layer methods, the 
problem to solve is essentially that of an inviscid flow 
whose boundary conditions are no longer known in advance. 
On the contrary, these become solutions of a '.system of 
ordinary differential equations on the boundary as, for 
example (3.7) (4.8) (4.9), relating the pressure p(xjO) 3 

the angular direction of the flow < ®(x, a ) as well as- the 



viscous thicknesses. After a few eliminations/ we may 
extract the relationship between yix.o) and p(x 3 o) actual- 
ly applied to the inviscid flow: 

4, e - At S ' , 4, (c 

x is measured along the boundary of the inviscid flow - 
Ai, A 2 j A 3 , 6 * depend on other differentail equations, on 
the initial values of the viscous thicknesses as well as 
on nonviscous variations & (x) and p(x). Equation (5.1) 
leads to highly nonlinear behaviors; Ai is reduced to zero 
during separations or reattachments, A 2 is reduced to zero 
at the critical points when these were not eliminated. The 
small scale term 6 * leads to phenomena of branching solu- 
tions, whose importance was reviewed and which lead to the 
inclusion of supplementary downstream boudnary condi tons. 


In the case where the integral boundary layer methods 
are replaced by local Prandtl viscous equations, the beha- 
vior of the viscous boundary, as a boundary condition of 
the inviscid flow, remains appreciably the same as (5.1), 
at least if we confine ourselves to an analysis of distur- 
bances around a solution 

( / K(*»] ■■ 

U a /•/£- %■] ( 


a relationship to compare with (4.2) (4.3) (4.11) (4.14). 

Use of the fundamental relationship (5.1) as a closing 
boundary condition of an inviscid flow calculation there- 
fore represents in a schematic form the coupling problem 
to be solved, in the case of thin viscous layers. Note 
that the viscous relationship (5.1) does not give the 
variables / 8(x,o)' and p(x,o) a symmetrical role; only the 
pressure has a function in the form of a derivative. This 
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major difference with the nonviscous problem undoubtedly 
is due to the thin layer aspect more than to the viscous 
approximations used. 

5.3. Methods Of Initial Supersonic Conditions 

The local function of relationship (5.1) initially 
led to an investigation of supersonic solutions to the 
strong coupling problem, in this case, the apparent simpli- 
fication of the initial conditions methods, proceding from 
upstream to downstream with simultaneous resolution of (5.1) 
was widely used by Crocco-Lees and numerous successors [32, 

33, 43, 59, 62 to 66] particularly in the case of a simple 
wave. However, we know that the unstable branching solutions 
result in an 'incorrectly defined problem, requiring a tech- 
nique of resolution by successive probes over an initial 
disturbance, so as to select, when convergence is reached, 
a solution having a behavior or stability which , is adapted 
to downstream conditions. The scheme in figure 16 shows 
that the initial disturbance e was conventionally produced 
on the pressure in subcritical cases, whereas it was pro- 
duced over the initial position of the pressure jump and 
thus enabling the solution to cross a critical point lo- 
cated farther downstream, in the most supercritical cases, 
and without singularities. 

In practice, such a resolution is extremely difficult 
and lacks generality. It is virtually unachievable as soon 
as-ithe region where the strong interaction begins, i.e. the 
disturbance region, is not known in advance, or as soon as /II 

successive strong interactions are present, a fortiori if 
they mutually interfere with each other. 
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5.4. Pseudo-Unsteady Methods 


The first important improvement was suggested by Werle- 
Vatsa 167 to 69] for simple wave type external supersonic 
flows in which the nonviscous calculation is reduced to a 
simple local algebraic relationship between the pressure 
and direction of the flow. By means of a slip coupling 
over the displacement surface, we obtain : 


p = - 6 5* 

r X XX 


(5.2) 


Werle-Vatsals basic idea was simply /to replace the coupling 
problem with a, fictitious problem of variation obtained : by re- 
placing the term 8* with [6* - 6*], and calculating the 

xx xx t 

asumptotic state when t -*■ » using a numerical method of 
alternating directions with two steps, the first being implicit 
in Cy 3 t ) , the second in (x 3 t). If we consider, for example, the 
equation of motion in x which, contains the only unsteady term: 
pu u x = py u y = + uy ^y (5.3) 


we obtain schematically to pass from time n to time (n+1 ) : 

I * e**p y w -I- ( -*”)] + [ £ */]y 

i 

I , , 

I [f“*X * r } 1 / [(K)" 1 - ( f**' - )] ]y 


(5.4) 


The first half time step leads to a conventional calcu- 
lation of the boundary layer where 8**, and therefore 
p is applied. It is resolved upstream to downstream, as 
the linearisation of the equations calls on the preceding 
time step introduced to avoid singular behaviors at the 
separations. The second half time step is reduced to: 


6 (8* ) n+1 - 5* ) n = h [ 5 * n+1 - 28* n+ % + 8 * n ] (5.5) 

— — At 


xx 


XX 


which leads to the simple resolution of the tridiagonal 
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n +2 

matrices for 6* , after discretization. The method 

was extended by Napolitano, Werle, Davis [69] to semi- 
unlimited incompressible external flows within the frame- 
work of laminar triple-deck equations. In this case, the 
albegraic relationship (5.2) should be replaced by a Cauchy 
integral (5.6) which makes the half-time step more complex 
(5.7): 


' f m Jx 

77 X-$ 


(cr. ± r 

77 J. m X. 


.zr 


rj 


X. % 


d > 


(5.6) 

(5.7) 


After discretization (5.7) leads for 6* n+ 1 to an algebraic 
system exhibiting a full matrix, but with diagonal domin- 
ance. A relaxation iteration is therefore necessary, in 
incompressible flows, to solve the second half time step 
by successive inversions of simple tridiagonal matrices. 


A quitd differerit pseudo-unsteady method was suggested 
by Briley, McDonald [70] and was used by Gleyzes, Cousteix, 
Bonnet [71] to calculate short transitional separation bulbs 
of leading edges, in incompressible flows, using a simpli- 
fied local strong coupling approach. The basic idea was 
to use real unsteady equations of the boundary layer, while 
retaining the steady inviscid flow equations. The resolu- 
tion is also pseudo-unsteady, to the extent Where the bound- 
ary layer is calculated from the nonviscous pressure field 
of the preceding instant, the inviscid flow being adjusted 
a posteriori as a steady flow under the most recently cal- 
culated displacement effect. This simple technique is 
still justified only by the obtainment of asymptotic steady 
solutions when £-»■<». 


5.5. Direct Or Inverse Methods Of Relaxation 

This principle of this class of methods lies on a more 
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rigorous formulation and on an extension of Conventional 
iterative methods, used to improve the weak coupling 
techniques at the cost of a few successive approximations 
over the displacement thickness. 

The strong coupling is obtained here for a rigorous 
convergence of the successive decoupled resolutions of the 
inviscid flow and unsteady boundary layers problems, while 
excluding any arbitrary smoothing technique. The stability 
of the iteration is controlled by sub-relaxation [3, 13, 39, 

72, 73] . The linkage of the segmented viscous and nonvis^ 
cous calculations isiquite convenient and leads to a broad 
generality of use in subsonic, transonic or supersonic states. 

The downstream influence on the upstream in supersonic flows /12 
may be obtained automatically [3, 13, 15], as in pseudo- 
unsteady techniques, by directly applying the appropriate 
boundary condtions downstream, whether strong successive 
mutually interfering interactions, or streams not reducing 
to a simple wave are involved. The counterpart of this 
generality appears, for certain cases, in the preliminary 
calculations required for the coupling iteration. However, 
experience has shown [3, 13, 15, 39, 77] that this coupling 
iteration may be easily compete with those required to cal- 
culate inviscid flows, in the case of subsonci or transonic 
relaxation methods. In this case, the coupling iteration 
appears as an explicit relaxation technique for the vis- 
cous boundary condition (5.1) applied to inviscid flows, 
and the increased computer costs associated with viscosity 
is limited. 

Considering that relationship (5.1) becomes identified 
with a boundary layer resolution, the coupling iteration, 
called direct or inverse, depending on whether we are sol- 
ving for the cycle (n+1) a direct or inverse type nonviscous 
problem: 
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Direct : 


(5.8) 


Inverse : 



Ay 

f Hi 

A 3 

S'Hi 


(5.9) 


For the direct method, if we assume that the inviscid 
flow is an operator which relates a pressure p n+ ^(x) to a 
given angular.vdistribution 9 % (x ) , the inclusion of (5.8) 
leads to operator 1 /" related to the iteration of the overall 
fixed point: 


9 v1 (x)m F[o'(x)] 


(5.10) 


Likewise, when A ^ may be reduced to zero owing to a 
separation, the resolution (5.8) should be avoided and, 
assuming that an inverse inviscid flow operator correlates 
| @'(x) with the data of p n (x), by adding relationship (5.9), 
integrated in x, we obtain an inverse fixed point iteration; 

p n+1 (x) = G [p n Cx)] (5.11) 

The adoption of one of the two techniques (5.10) (5.11), 

or even their joint use in an alternating' combined method 
(5.10) or (5.11) as a function of the predefined regions 
on boundary Ox, were used by Lock [39], Lineberg, Steger 
[74], Le Balleur [61, 72], Carter [73], Melnik, Chow, Mead 
[26, 29], among many other authors. In all of these cases, 
the method was stabilized by a large sub-relaxation u, in- 
dependent from Xj determined without any other criterion 
than the previous numerical experience: 

j e'"(z) . 9'(x) * « {F[e •(*)] - o'(*) j (5.12) 

j 

The lack of generality of the previous approach was, 
however, eliminated by the present author through a- study 
of approximated stability [13] permitting a local calculation 
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of Ui(x). The study shows thht operator F of (5.10), 
linearized around e"(x.) . , appreciably implies that the 
harmonic distributions of disturbance [q ( X ) - &“l x ) J are 
specific functions and that the corresponding approxi- 
mation of the specific complex values of F is yo, that 
of G being yi with: 


/■o ■/*/ * 1 


ccBS* 

fiTti 


(5.13) 


M designates the local Mach number of the nonviscous flow, 
a the harmonic frequency, B is a parameter of viscous 
form related with B' in (4.11), or (A 2 /A\) in (5.1). 

The. highest specific value yo (x) or y i(x) makes it pos- 
sible to determine u(x)i hiqhly deoendend upon x like 
yo or yi (see [13]), while a local over-relaxation is not 
te be exluded. 


If i?( a) and: / (*) designate the real and imaginary parts 
of y 0 for example, the stability of (5.12) will be obtained 
if: 


to (x) 4 2 (* mjx ) 


Vft - 


1.R 

r 



(5.14) 


In the direct case (5.10), a discretization step, 
small Ax increases the highest frequency a of the possible 
disturbances and reduces the value of u(x) in the direct 
problem (5.10), while to must be zero at separation or re- 
attachment points where A 2 = 0 3 B ->«>, y 0 In inverse 

mode (5.11), however, the cancellation of yi leads to these 
points at (o = 1. 

Numerically, the local calculation of w(x) shows that 
the direct mode relaxation (5.10) is well conditioned for 
attached boundary layers, the inverse relaxation on the 
pressure (5.11) is also near separated layers or when it 
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is close to being near separated layers. 


5.6. Implicit Linearized Method 


Besides the dual nature of the preceding behavior, 
physically related to a dominant influence of the inviscid 
flow or the boundary layer for the pressure deterimination, 
depending on whether the boundary layer is more attached 
or more separated, the fact that there is a stronger ^re- 
course to a sub-relaxation when hx is smaller shows the 
explicit nature of the" method (5.10) as well as the pro- 
bable advantage of a more implicit method. The idea of 
a Newton method for solving the discrdtized equation (5.10) 
was suggested by Burne, Rubbert, Nark [60] . Accordingly, 
for each iteration, we must determine the inviscid flow 

71 71 

and the boundary layer influence matrices. They 

are expressed if we select for variables p(x) and S*(x), 
discretized in p. and 6* in the inviscid flow, in p. and 
<$| in the boundary layer: 


l T- % (*-%) 


(5.15) /13 


For iteration (n+1) we may therefore adopt the solution 
of the linearized coupled problem which is expressed in the 
form of a matrix: 


I f” = „ /* * (i.p^j' i [p x (r n . i**) , 

r", r\ 3'[(r-r)>(r'-r)] 


(5.16) 


If the nonlinear effects on P n - • and B. . are weak, 

" 0 I'd 

the convergence must be rapid. Each iteration is still 
complex, requiring the inversion of a full matrix [I - P n B n ] 
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and above a 1*1 ah estimate of the influence matrices P n and 
B n , a priori a costly evaluation within the framework of 
finite differences methods, particularly in transonic flows. 
Little information is currently available about this method, 
although its use is indicated by Thiede [54] . A very sim- 
ilar technique is described by Arieli, Murhpy [76] for 
pseudo-direct calculations of the boundary layer. 

5.7. Semi-Inverse Relaxation Methods 

The the practical standpoint, the complexity of the 
method described above as well as the mixed type relaxation 
iterative methods (.direct or inverse by regions) leads us 
to suggest [3, 13 to 15, 75] a new relaxation iteration for 
separated boundary layer regions or those on the verge of 
reaching this state. This method may be qualified as semi- 
inverse in that only the boundary layer problem is processed 
in inverse mode, the given data being either the direction 
of the external f low &'(x) ,j or the displacement thickness 
<$* Cx) . As shown in figure 17, &*(x) gives a double pres- 
sure predictor for each iteration, one given by the inviscid 
flow p Cx) , the other given by the boundary llayer p(x) , the 
problem. being to correct iteratively so that p n (x) 

j will converge near pCx). 

An approximated stability analysis of the direct or 
inverse coupling iterations described in paragraph 5.4 
enabled us to offer a solution to this problem. We may 
notice that if the correction selected - &'(x) J is such 

that the disturbance [p n+ 1 (x) - p™Cx)] in the inviscid flow 
is exactly IpCx) - p n Cx)], the new distribution . is 

identical to the result of an inverse calculation of the 
inviscid flow subjected to pCx). The chainage of such an 
inverse calculation simulated for the inviscid flow and a 
true inverse calculation for the boundary layer, therefore 
reproduces iteration (5.11) and may be stabilized like it. 
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by a local calculation of a relaxation coefficient w(x) , 
according to (5.13) (5.14). The analysis must still be 

completed by an appropriate technique in order to solve 
the problem of small inverse disturbances, of the Prandtl- 
Glauert type, enabling JS'%)- &( X )J to be calculated from 
[ p (x ) - p n (x )] . We have shown [13] that a harmonic anal- 
ysis makes it possible to reduce each frequency mode a 
to a local algebraic relationships whose association with 
with the relaxation (5.13) (5.14) leads us to the following 

corrections, if U n (x) and U(x) designate the velocities re- 
lated to the pressures p n (x) and p(x): 

Subsonic points: tfBfTrF f f dU_ j_ do '* 1 

0 -6 . T & U' dx J (5.17) 

Supersonic pbints $*• - -V ^~rr / 

* 0 - m 0 / T dz‘ U % dz‘ J 

by adopting a max a = ( ir/A ) , the dependency of the preceding 
local corrections on 8 X with respect to a are finally elim- 
inated. Among other advantages, the semi-inverse iteration 
(5.171 requires only a detection of the coupling errors on 
the pressure gradient. As a result, it can easily be alter- 
nated with the direct iteration (5.12) as a function of the 
local form parameter of the viscous layer, without ever lead- 
ing to discontinuous pressure distributions during the course 
of the convergence, for any number of regions resolved in in- 
verse mode, or for any modifications made during the itera- 
tions . 

For example, this combined method was used to calculate 
the supersonic separation on a compression bar, in a laminar 
or turbulent flow, on the assumption of a simple wave, by 
solving a well-defined problem directly with a zero pressure 
gradient at the last calculation point, figures 18 and 19. 

It was also applied to calculate symmetrical transonic air- 
foils in theories of. small disturbances with shook or lead- 
ing edge separation and to capture numerically the strong 
interaction at the base of a shock over a fine mesh, figures 
20 to 22. 
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A semi-inverse slightly different correction was 
suggested by Carter [78] more recently. For each itera- 
tion, the correction is made of the displacement thick- 
ness 6* Cx) and not of the angular direction of the flow. 
In the incompressible case, for example, the correction 
is expressed in a local algebraic form: 


S* n+1 = ft 8* n [ jjn ] (5.18) 

in which !1 is a relaxation coefficent, even over-relaxa- 
tion coefficient (1 < $2 < 2) . This method was recently 
used by Kwon, Pletcher [79] , Whitfield et al. [80] . 


It is interesting to note that there is an analogy /14 

between (.5.18) and (5.17) in the subsonic case. If we 
assume, for example, that the coupling is obtained by 
inclusion of the displacement thickness, or & » (d&*/dx ) , 
the derivation of (5.18) directly gives: 


_ Jill » nt* _L 1 

dx dx L u dx- U* dx J 


(5.19) 


The rigorous identity of the formulations (5.17) (5.18) 

should therefore be obtained for the incompressible case, 
if ft = [Bkx/-n§*B-hx)] . The analogy (5.17) (5.19) also 

shows the ambiguity of the notion of over-relaxation or 
subrelaxation in (5.17) (5.18), the subrelaxation in 
figure 17 occurring in (5.17) only at the concept level, 
the over-relaxation ft of (5.18) being weighted by the 
multiplicative factor 5*. 


Finally, let us recall that the use of non-automatized 
semi-inverse type techniques was first mentioned by Kuhn 
[81, 82]. 
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5.8. Serai-implicit Methods 


Recent trials have been carried out to develop methods 
which are partially free of the explicit nature of the re- 
laxation techniques described in paragraphs 5.5 and 5.7, but 
leading to simpler solutions than Newton's method in para- 
graph 5.6. 


A 1 first analysis was suggested by Veldman [83] in the 
case of unlimited incompressible flows for which the rate 
of disturbance induced on the viscous layer is given by a 
homologous Cauchy integral of (5.6), if the coupling is ob- 
tained by a distribution of sources on the wall: 

(5.20) 

J-*> z a y 

After the discretization, the velocity distribution 

u. or pressure distribution p. becomes associated, by virtue 

of (5.20) with the displacement thickness distribution 

by a linear form whose coefficients A. . , pr edontinantly dia- 

^ <7 

gonal, can be calculated directly. We may express: 

Inviscid flow: p n . +1 - A . . &**+* ~ E A. . 5 * n+1 + X (5.21) 

* *3 * j<i *3 3 j>i ! ' '■ 


crete 
a tor 


Schematically, the boundary layer appears in a dis- 
form like an appreciably lower ^triangualr matrix oper- 



Boundary Layer: 


st n+1 


r.n+1 n+1 v u n+l 

B . . p . = .4 . B . 


n+1 


+ .Y. . 

J 




(5.22) 


It is clear that the iteration indices n and (n+1) 
were selected in (5.21) and (5.22) so that a resolution 
worked: out from upstream to downstream may implicitly 
solve the first members in station i, the second members 
being known, according to a Gauss-Seidel type technique. 
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Like the Gauss-Seidel technique/ the method (5.21) (5.22) can 

also be over-relaxed. 

; * • i 

A second semi-implicit analysis developed by Wai, Yoshihara 

[84] for calculations of transonic airfoils in a small distur- 
bance approximation. The potential problem: 

' (I- M* . ><p x )fr c * - 0 (5.23) 

; • £(*)« (f&x) (5.24) 

is solved, in the inviscid flow where (x) by the Murman-Cole 
relaxation method/ implicit column-wise along y, semi-implicit or 
implicit in x, depending on the decentering applied t ° <P X and <£ .' 

In the inviscid flow, the boundary layer behavioral relationship 
(5 . 1) (jives: 

i 

4> ,,(x» 0) = fix) + b(x) <p i c(x) (5.25) 

y xx 

a relationship in which b(x) and o(x) should be calculated by the 
other boundary layer equations. If this calculation is performed 
from upstream to downstream while scanning the columns at the same 
time for <j>, it is clear that (5.25) (5.23) may be solved simultan- 

eously in serai-implicit mode, exactly like for inviscid flows, with 

the same discretization for & and <]> .T The main advantage is that 

y xx 

the interactive treatment (5.25), even when discreticized in semi- 
implicit mode for <P f the improvement being that of a a Gauss-Seid- 
el technique, in analogy with (5.21) (5.22). 

The simultaneous relaxation of the integral boundary layer 
equations and of the stream function, through successive implicit 
scannings by columns, was also indicated by Moses, Jones, O'Biien 

[85] in the case of separated subsonic flows. It is used by Ghose, 
Kline. [86] in the calculation of diffusers. 

6. SYNTHESIS EXAMPLE - CALCULATION OF AIRFOILS 

Among the methods of calculating compressible flows of viscous 
fluids past airfoils, using coupling methods like those used in 
the approximations in paragraphs 4 and 5, 
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we may mention those of Bavitz [37], Bauer, Korn [38], 

Lock, Collyer, Firmin, Jones [39] , Melnik, Chow, Mead 

[29], Wai, Yoshihara [84], Nandanam, Stanewsky, Inger /15 

[40], Le Balleur [15], such a list not being limitative. 

We are not interested here in the comparative performan- 
ces of the methods from a practical standpoint, most still 
being the subject matter of current developments and gen- 
erally providing more accurate results than currently us- 
ablfe Navier-Stokes solvers [87] . What we are looking for 
here is simply an example of a synthesis calculation, dir- 
ectly usable in applications, combining various viscous 
interactions mutally influencing each other (trailing edge, 
shock-boundary layer interaction, separation) , and reveal- 
ing the development of methods for an effective treatment 
of strong couplings and those offering a complete enough 
methodology to completely describe the viscous flow at 
high Reynolds number, in a direct extension of numerical 
techniques developed for inviscid flows. 

6.1. Bavitz, de Bauer, Korn et al. Methods 

In these first methods of calculating airfoils in 
viscous flows, the approximations used for the coupling are 
still close to the weak coupling methods. The wake is still 
not taken into account. The displacemement effect of the 
boundary layers is obtained by altering the airfoil design 
geometry after adding the displacement thickness. The 
stabilization and convergence of the coupling iterations 
are relatively uncertain and incorporate smoothings or ex- 
trapolations of 8*(x), an empirical processing of the trail- 
ing edge region being used as a compensation in the Bavitz 
method. The inviscid flow is calculated using the poten- 
tial Garabedian, Korn method. The nonconservative form is 
mos often preferred owing to a certain error compensation 
bectaus^ the viscous effects at the base of the shock are 
not taken into account, the mesh being selected so as to 
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extend the compression under the shock far enough to avoid 
any separation in the boundary layer calculation [37] . 

6.2. Collyer, Lock et al. Methods 

The investigation of more complete methods for strong 
interacting phenomena was initiated by Lock and his asso-’ 
ciates by first solving the inviscid flow using small dis- 
turbances transonic techniques/ figure 23, then using Gara- 
bedianls land Korn's potential method, and this was inde- 
pendently of the practical performances of the analyses men- 
tioned above. Figure 23 recalls that the wake effects were 
introduced in the form of a normal velocity jump, corres- 
ponding to the displacement effect (3.5) and in the form of 
a tangential velocity or pressure jump associated with the 
wake curve effect of the type (4.13). The curve of the 
stream lines under consideration is still a mean curva- 
ture. Furthermore, in the case of the complete potential 
equation, the wake conditions are applied to a mesh line 
and not to the viscous wake itself. The separations are 
not processed. 

Compared to the preceding methods, the Imain improvements 
besides treating the wake are that the normal velocity on 
the airfoil wall is processed to determine the displace- 
ment effect without altering the design geometry, the sim- 
ultaneous iteration on the viscous effects and on the 
potential calculation, the replacement of the smoothing 
techniques on 6* with a uniform viscous sub-relaxation, 
determined empirically without altering the solution when 
convergence is reached. The latter improvement consisted 
of simulating as close as possible the pressure jumps of 
the shock waves provided by the potential calculation, by-, 
weighting the conservative and nonconservative numerical 
techniques. No special viscous processing was introduced 
for the shock-boundary layer interaction. 
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6.3. Melnik, Chow, Mead Method 


Like the method just described, this anlysis is 
applied to airfoils free of separations, based on the 
Green integral method ofr boundary layer and wake cal- 
culations, and does not introduce any special processing 
for the strong interaction at the base of the shock waves. 
Figure 24 recalls that, excluding the trailing edge region, 
the viscous boundary conditions applied to the inviscid 
flow are similar to the method of Lock et al., the condi- 
tions on the wake always being applied to a slightly dif- 
ferent line than that of the real wake. 

The improvement of the methods consists of using the 
Jameson potential calculation in a conservative form, and 
using a sophisticated processing of the strong turbulent 
interaction at the leading edge: the latter is based on an 
asymptotic multi-deck analysis, figure 25, and guarantees 
a rational behavior of the radoel at the infinite limit of 
the Reynolds number. From the practical standpoint, we 
should notice in figure 25 that the viscous calculation 
model incorporates a non-irrotational normal pressure gra- 
dient in the external layer. This is expressed in the po- 
tential calculation by a local modulation of the tradition- 
al coupling boundary conditions at the leading edge, figure 
24. it generally leads to overvalued pressure jumps on the 
wake in the immediate vicinity of the trailing edge. 

6.4. Nandaraan, Stanewsky, Inger Method. 

This method is simpler thant the two preceding ones, 
except for the processing of the shock-boundary layer 
interaction. It uses the Jameson potential method, in a 
conservative and nonconservative form, as well as an inte- 
gral method for the boundary layer. The wake is not computed. 
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The separations are not processed. The coupling is obtained 
by altering the design geometry of the airfoil by adding 
a displacement thickness subjected to a large smoothing, 
particularly in the shock region. 

The method would therefore be comparable to those 
described in paragraph 6.1., if a more sophisticated proces- 
sing of the shock-boundary layer interaction were not incor- 
porated. The latter appears as a separate module whose role 
is to determine as accurately as possible the state of the 
boundary layer after the shock, and to provide a more real- 
istic estimate of the pressure on the wall in the shock re- 
gion. This module is based on Inger's non-asymptotic multi- 
deck analysis [88] . The results are better if the compres- 
sion applied to the Inger module, theoretically that of a 
straight shock, is that of a maximum deflecting oblique /16 

shock. 

6.5. Wai, Yoshihara Method 

The method is a recent development and is used for the 
potential flow of a simple small trahsonic small disturbances 
calculation. The boundary lyaer and the wake are calculated 
using an integral method, including separations. The first 
feature consists of a simultaneous resolution of the poten- 
tial small disturbances equations and of the viscous bound- 
ary condition C5.1), during the relaxation by column' asSa- I 
ciated with the Murman-Cole technique. This process, eli- 
cited in paragraph 5.8, is semi-implicit in that the scan- 
nings per columns must be repeated to assure the convergence. 
Conversely, the interactive processing of the viscous calcu- 
lation eliminates the singularity problems at the separation. ' r '~ 

The second feature of the method is to replace the 
modelization, called viscous ramp modelization, with an 
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integral calculation of the boundary layer in the shock 
region, in the case of separation, in the absence of re- 
solving the strong coupling with a fine enough mesh. The 
modelization schematically consists of establishing a ramp 
for the displacement surface with an angle which corres- 
ponds to a predetermined recompression level on a shock 
polar, where the empiricism leads to the selection of an 
intermediary value between the sonic deviation and the max- 
imum deviation [89] . This modelization of 6* is necessary 
not only for improving the comparison with the experiment, 
but also to assure the convergence of the calculation. 

6.6. Le Balleur Method 

The method uses for the inviscid flow the potential 
calculation of Chattot, Coulombeix, Tora§ [90] which solves 
the equation in a conservative form, with the Jameson arti- 
fical viscosity. The boundary layer, laminar or turbulent, 
as well as the wake are calculated with an integral method 
including the reverse flows, the wake being calculated with 
or without the dissymetry effect. The matching formulation 
of paragrph 4.3 for the inviscid flow leads to boundary 
conditions on the wall and on the wake which are similar to 
those of Lock et al. or of Melnik et al. . Yet the geometry 
of the wake line is here still periodically dapated during 
the calculation to coincide with the mean viscous wake line. 
The method does not provide any special processing of the 
leading edge, but strictly observes here the strong coupling 
assicated with the matching formulation. The curvature ef- 
fect applied (4.13) incorporates the averaged estimates of 
K(x), which are distinct for the upper and lower half-wakes. 

The coupling is achieved by a direct' or semi-inverse 
relaxation iteration, as a function of the instantaneous 
local form parameter of the viscous layer. This iteration, 
described in paragraphs 5.5 and 5.7, is comparable to that 
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used for calculating the inviscid flow, figure 26. The 
pressure gradient along the wall is estimated in the in- 
viscid flow using centered schemes for subsonic flows, 
and decentered downtream for supersonic flows. The bound- 
ary layer transition at the wake is assimilated with a con- 
tinuous process for the inviscid flow as well as for &*(x) 3 
but discontinuous for the viscous calculation. The method 
automatically captures the modification brought to the 
pressure distributions by the. viscosity in the immediate 
vicinity of the trailing edge, provided that the meshes 
are approximately- -the local thickness of the viscous layer, 
whether there is a separation or not. 

Figures 27 to 32 show the use of the method at low 
velocities and at high incidences, for the NACA 0012 airfoil. 
On this airfoil, the calculation approaches the stall inci- 
dence, with a 70%r/upper -surface ^separation from the chord-, 
and a reattachment in the wake. In the present state of 
of the method's development, this calculation at an inci- 
dence of 16° is possible only by observing the wake dis- 
syraetry, figures 27 and 31 show the aspect ratio .change 
associated with its mean line. Figures 39 to 41 show by the 
plotting of the stream lines and of the iso-Mach lines that 
the calculation method restores not only the distant invis- 
cid flow, but also the mean viscous flow in detail. 

Figures 33 and 34 show the use of the method in sub- 
critical transonic cases on the RAE 2822 airfoil, tested 
by Cook, McDonald, Firmin [91J . In supercritical- transonic 
cases, in contrast to the calculations made in the small 
disturbances approximation, figures 20-21, it remains cur- 
rently impossible to numerically capture the strong inter- 
action at the base of the shock waves on a fine mesh adap- 
ted to the local boundary layer thicknesses, owing to the 
riiimerical difficulties of solving the potential equation 


46 



on such meshes. Without this complete resolution of the 
local viscous coupling, the compression under the shock 
is overestimated by the conservative method applied to 
the inviscid flow. 

To correct this problem, an approached method was 
defined for the shock-boundary layer interaction. As fig- 
ure 36 schematizes, this method simply consists of solving 
a local fictitious interaction on an enlarged scale e 2 , 
compatible with the mesh. This is a homothetic image of 
the real scale interaction ei, obtained by simple expan- 
sion of the local thickness 6 of the viscous layer in the 
ratio (e 2 /si), yet without altering the Reynolds number 
ES . The strong coupling made over the fictitious inter- 
action region e 2 leads, on the one hand, to a more realis- 
tic interactive estimate of the boundary layer thickness, 
and on the other hand, to an expected reduction in the re- 
compression level, as shown in figure 37. In contrast, the 
processing is essentially an approximation in that the al- 
teration undergone by the inviscid flow in region e 2 is 
inconsistent with the overall calculation scale C. For 
this, a quality numerical solution will therefore be j 
maintained only if c 2 << and even e 2 •£< L 3 L represent- 
ing the extent of i: the supersonic pocket. It may be observed 
in figure 38 that, thanks to i thee. single viscous coupling 
effect, the calculation predicts a Mach number, after a 
shock which is slightly less than unity , which is closer 
to the experiment. It may also be observed that the numer- 
ical structure of the shock of the potential calculation is 
replaced by a continous viscous compression, which is rela- 
tively spread out owing to the mesh used, figure 37. Fig- 
ures 42 to 44 show the appearance of the iso-Mach lines and 
of the stream lines in the calculated visedus fldwi. 
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7. CONCLUSIONS 


/17 


A spectacular effort has been made over the past few 
years in the inviscid flow region with a view to reinforc- 
ing the effectiveness and generality of the calculation 
techniques. It is important to stress first the major 
valorization resulting' :in practice from the coupling of 
these techniques with complementary viscous calculations, 
provided, however, that too restrictive approximations of 
the weak coupling are avoided. Changes in the methods of. 
predicting local and overall wing and airfoil character- 
istics shows this trend. 

In this light, we may note the conjunction of progress 
made "Of! i.an analytical natute on interactive viscous calcu- 
lations of fundamental phenomena such as trailing edge flows, 
thin layer separations, boundary layer - shock interactions, 
with progress made of«>.a numerical .mature, consisting mainly 
of developing new original strong coupling techniques. These 
improvements have made it possible to go beyond the stage of 
processing specific or local phenomena, and to synthesize 
computer means adapted to the complex flows encountered in 
industrial applications which compete with overall methods 
for direct resolutions of the viscous flow. 

Yet, it is clear that the two approaches to phenomena 
of viscous interaction, the overall approach and the strong 
coupling approach, are actually more complementary than com- 
petitive. The current selection of relatively simple vis- 
cous approximations in the coupling methods is actually 
guided only by the concern to use coarser discretization 
meshes than those which are indispensable in the overall 
approach for numerically solving viscous layers and sub- 
layers. Coupling methods using simple viscous equations, 
and particularly those of Navier-Stokes type calculations. 
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produce in practice the minimum mesh densities compatible 
with the numerical resolution of strong coupling phenomena. 
Such densities, are between those used in inviscid flows 
and those of the Navier-Stokes type calculations, and lead 
to an economy which may be taken advantage of for calcula- 
ting more complex flows, incorporating for example multiple 
viscous interactions or three-dimensional phenomena. 

For less restrictive meshes, the progress made in the 
overall approach has obvious effects on the calculation of 
viscous regions in the method of coupling by regions. Fin- 
ally, it is also conceivable to use Navier-Stokes solvers 
in the :;<i6up ling methods replacing current inviscid flow cal- 
culations and on analogous meshes. Such a technique would 
make it possible, in three-dimensional flows, for example, 
to capture the macroscopic viscous phenomena by a local re- 
solution of the Navier-Stokes equations with a turbulence 
model, whereas the concentrated phenomena of thin layers 
near the walls, their interaction with the external flow 
as well as their impact on the generation of vortex sheets 
would continue to benefit from legfe l.costly processings, 
through extensions of integral interactive boundary layer 
methods currently developed for airfoils and wings with a 
large aspect ratio. 

From the specific standpoint, then, a few salient ideas 
or requirements may be extracted from current coupling tech- 
niques: 

(11 the implementation of strong interaction coupled 
calculations presents two relatively different methods. The 
first uses a concept of regions linked by their boundary 
conditions, with the possibility of directly incorporating 
existing modules of the Navier-Stokes type. The second is 
based on the notion of overlapping viscous and nonviscous 
calculation regions, with the possibility of developing 
original methods for solving Navier-Stokes equations using 
an equation scattering technique. 
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(ii) The necessary replacement of weak coupling 
approximations with a strqng' rigorous coupling leads to 
an imperative requirement to develop specialized numer- 
cical methods/ raising problems that are at least as com- 
plex as those associated with numerical techniques for 
inviscid and viscous flows. This requirement became ob- 
vious for trailing edge, separation and shock-boundary 
layer problems. The analysis is especially not very de- 
veloped for a coupling between the Euler and Navier-Stokes 
equations . 

(iii) The processing of boundary layer - wave shock 
interactions seems to be one of the most delicate problems 
of the calculation techniques, not only due to the error 
range attached to the turbulence mddelizations , but also 
due to the reduced scales which should be used locally, 
even for calculating inviscid flows in the case of coupling 
method using miiiiimal meshes. Similar problems of multiple 
resolution scales are presented in viscous layer separations 
in the vicinity of leading edges. 

. (iv) Among the improvements made in calculating air- 
foils, the replacement of the Euler equations for the po- 
tential equation shall be examined, particularly for higher 
transonic flows, in the event of a separation generating 
under the shock, with the ultimate benefit of a direct ap- 
plication to internal flows. Calculation-experiment com- 
parisons brought to light from the corrections of error 
ranges on windtunnel walls, should also be highly desirable 
for refining existing methods of calculation. 

.. (.vi Calculations of airfoils in extreme flows, near 
stall or separation states under the shock, for exmaple, 
reveal the very strong sfensitivity of the overall flow 
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to the precision of the viscous layer calculation from 
the very origin of the latter. The development of tur- 
bulence models is therefore very important, even in the 
simplified forms encountered in integral methods. 
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Figure 1 - Computation Methods In Viscous Fluids. 


60 




Fig 4 - Overlapping 
Coupling Approach 


Key: a -Methods; b-Over lapping meshes; 
c-Strong coupling over ft. 
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Figure 5 - Asymptotic Weak Coupling 
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Fig 8 - Matching formulation: Scheme of vis- 
cous and non-viscous solutions for a 
given x-axis. 
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Fig 9 - Laminar supersonic ramp - 
Asymptotic behavior of the 
interacting Prandtl equations, 
according to Burggraf, Rizetta, 
Werle, Vatse, ref. [27], 
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Fig 10 - Critical points in the 
Crocco-Lees sense. 

Key: c-Uns table branchings; 
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Fig 11 - Supercritical - Subcritical 
jumps. 
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Fig 12 - Supersonic separation over a ramjp. Diagram 
of the branching solutions for ain initially 
subcritical layer. 



Fig 13 - Supersonic separation with the coupling model 
for adjoining regions: a) subcritical; b) su- 
percritical. 



Fig 14 - Coupling of overlapping regions:' nonviscous 
schematization of a supersonic turbulent 
separation. Pressure distribution at the 
wall. 
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Fig 15 - Numerical instabilities of 
conventional couplina iterations on 

a flat plate in an unsteady flow as a function of Ax & At. 



effect 


Fig 16 - Methods of initial conditions in supersonic flows. 



Fig 17 - Numerical method for a semi-inverse coupling 65 
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• Experiment: Lewis-Kubota-Laes 

N.S. calculation: Carter 
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Present coupling method 


Fig 18 - Turbulent supersonic 
separation on a ramp. 


Fig 19 - Laminar supersonic. 

separation on a ramp. 



Fig 20 - Strong coupling calculation or a turbulent transonic 
separation under a shock of a circular symmetrical 
airfoil of 18%. M = 0.788 - R = 4 x 10. 
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Fig 21 - Shock-boundary layer interaction region 

in a strong coupling calculation in small 
transonic disturbances. Viscous thicknesses 
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Fig 22 - Symmetrical circular airfoil of 18% 
in small disturbance theory. 
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Fig 23 - Coupling calculation of airfoils 

in viscous flows: boundary condition 
of the inviscid flow, in small dis- 
turbances, according to ref. [77]. 



Fig 24 - Coupling calculation of airfoils in 
viscous flows: boundary condition 
of the inviscid flow according to 
Meinik, Chow, Mead, ref. [29]. 
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Fig 25 - Nonseparated turbulent trailing 

edge: asymptotic model by Melnik, 
Chow according to reference [29] . 
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Fig 26 - Mixed relaxationumethod : direct or 
semi-inverse . 
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Fig 27 - Mesh with positioning of the wake in 

(M = 0,776, a = 75° ,R = 2x Uf). 


the symmetrical 
approximation 
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Fig 30 - Viscous effects at the trailing 
edqe and pressure plateau. Coupled turbulent 
calculation, symmetrical wake. 
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(M = 0.116, a = 76°, R = 2x 10*). 

Fig 31 - Mesh with wake positioning in the 
dissymetrical approximation. 
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Fig 32 - Calculation of the 
NACA 0012 at high incidence 
with dissymetrical wake. 


Fia .33 - Calculation of the RAE 
a i r foil in a subcritical flow. 
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Fig 35 - Mesh in a supercritical transonic flow with 

- symmetrical wake. 
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Fig 36 - Approximate calculation of the shock - 
boundary layer interaction if 
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Fig 37 - Effect of the shock/boundary layer 
interaction. 


*Experiment 



2822 


Fig 38 - Calculation of the RAE 2872 in a 
supercritical flow. 
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Fig 42 


Profit RA E 2822 supercrttlque (M = 0, 732, a = 3°, R = 6,5 x 1(f). 

Isomach lines\ in an inviscid flow — RAE supercritical 
^ | 2822 airfoil 

X N \ ! 

/ ■ \ \ 



/ Fig 43 - 


Isomach lines in a viscous flow - supercritical 

' „ RAE ,2822 airfoil 

Profit RAE 2822 supercrttlque (M =0,732, a = 3 , R =6,5 x 10°). 


Profit RAE 2822 (M = 0, 7 32, a = 3°,R = 6, 5 x 10°). 

Fig 44 - Streamlines in an inviscid flow 





